Smoothing and Change Point Detection for Gamma Ray Count Data   总被引：1，自引：0，他引：1  

Daniela Leonte  David J. Nott  William T. M. Dunsmuir 《Mathematical Geology》2003,35(2):175-194

Gamma ray detectors are used to measure the natural radioactivity of rocks. For a number of boreholes drilled at a site the gamma ray detector is lowered into each borehole and counts of gamma ray emissions at different depths are recorded as the instrument is gradually raised to ground level. The profile of gamma counts can be informative about the geology at each location. The raw count data are highly variable, and in this paper we describe the use of adaptive smoothing techniques and change point models in order to identify changes in the geology based on the gamma logs. We formulate all our models for the data in the framework of the class of generalized linear models, and describe computational methods for Bayesian inference and model selection for generalized linear models that improve on existing techniques. Application is made to gamma ray data from the Castelreagh Waste Management Centre which served as a hazardous waste disposal facility for the Sydney region between March 1974 and August 1998. Understanding the geological structure of this site is important for further modelling the transport of pollutants beneath the waste disposal area.  相似文献   

TEM study of a special grain boundary in a synthetic K-feldspar bicrystal: Manebach Twin     

S. Heinemann  R. Wirth  G. Dresen 《Physics and Chemistry of Minerals》2003,30(3):125-130

 A monoclinic KAlSi₃O₈ feldspar Manebach twin boundary was synthesized by diffusion bonding and examined using high-resolution transmission electron microscopy. The sharp (001) twin boundary is straight and free of strain. The boundary width is smaller than d₀₀₁. There is no rigid body shift observed at the twin boundary, and the feldspar structure is arranged symmetrically across (001). The twin boundary structure consists of bridged tetrahedral crankshafts, which are characteristic of the feldspar lattice. The grain boundary structure is in good agreement with the geometrical model of Taylor et al. (1934). The grain boundary composition of K_1/2H_1/2AlSi₃O₈ differs from their model. Received: 13 February 2002 / Accepted: 24 December 2002 Acknowledgements We thank M. Rühle, S. Hutt, J. Mayer, A. Strecker and U. Salzberger at MPI, Stuttgart, for their support and valuable advice in preparing TEM sections of bicrystals.  相似文献   

Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite     

E. J. Palin  M. T. Dove  S. A. T. Redfern  C. I. Saniz-Díaz  W. T. Lee 《Physics and Chemistry of Minerals》2003,30(5):293-304

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.  相似文献   

耦合海-气环流模式中双热带辐合带现象及其热收支分析   总被引：7，自引：2，他引：7       下载免费PDF全文

李江龙  张学洪  俞永强  戴福山 《气象学报》2003,61(1):39-51

文中研究了耦合海-气环流模式中的双热带辐合带(Double ITCZ)现象,并对这一现象的成因从海洋热量收支的角度进行了初步分析。Double ITCZ现象是在热带太平洋赤道南北两侧各出现一条ITCZ的现象,这是耦合海-气环流模式中的较为普遍的一种异常现象,与实际气候中出现的Double ITCZ现象并非指同一问题。文中对比观测和模式结果,指出了Double ITCZ现象的主要特征,针对它的出现过程进行细致分析,再利用模式输出的热量收支各项进行统计,得到了从海洋热收支分析得到的海表温度变化原因。与观测到的正常模态相比,Double ITCZ是一个异常的模态,它的特征突出地表现为西太平洋暖池区的降温和东南太平洋10°S附近的升温。海洋热量收支分析表明,暖池区的降温主要是受到扩散的作用,而表层平流的异常增暖在决定异常辐合带区升温过程中占第一位的作用。  相似文献   

The identifiability of parameters in a water quality model of the Biebrza River, Poland   总被引：1，自引：0，他引：1  

Marcel van der Perk  Marc F. P. Bierkens 《Journal of Hydrology》1997,200(1-4):307-322

The identifiability of model parameters of a steady state water quality model of the Biebrza River and the resulting variation in model results was examined by applying the Monte Carlo method which combines calibration, identifiability analysis, uncertainty analysis, and sensitivity analysis. The water quality model simulates the steady state concentration profiles of chloride, phosphate, ammonium, and nitrate as a function of distance along a river. The water quality model with the best combination of parameter values simulates the observed concentrations very well. However, the range of possible modelled concentrations obtained for other more or less equally eligible combinations of parameter values is rather wide. This range in model outcomes reflects possible errors in the model parameters. Discrepancies between the range in model outcomes and the validation data set are only caused by errors in model structure, or (measurement) errors in boundary conditions or input variables. In this sense the validation procedure is a test of model capability, where the effects of calibration errors are filtered out. It is concluded that, despite some slight deviations between model outcome and observations, the model is successful in simulating the spatial pattern of nutrient concentrations in the Biebrza River.  相似文献   

Multivariate nonparametric resampling scheme for generation of daily weather variables     

B. Rajagopalan  U. Lall  D. G. Tarboton  D. S. Bowles 《Stochastic Environmental Research and Risk Assessment (SERRA)》1997,11(1):65-93

A nonparametric resampling technique for generating daily weather variables at a site is presented. The method samples the original data with replacement while smoothing the empirical conditional distribution function. The technique can be thought of as a smoothed conditional Bootstrap and is equivalent to simulation from a kernel density estimate of the multivariate conditional probability density function. This improves on the classical Bootstrap technique by generating values that have not occurred exactly in the original sample and by alleviating the reproduction of fine spurious details in the data. Precipitation is generated from the nonparametric wet/dry spell model as described in Lall et al. [1995]. A vector of other variables (solar radiation, maximum temperature, minimum temperature, average dew point temperature, and average wind speed) is then simulated by conditioning on the vector of these variables on the preceding day and the precipitation amount on the day of interest. An application of the resampling scheme with 30 years of daily weather data at Salt Lake City, Utah, USA, is provided.  相似文献   

Variation of the plasma turbulence in the central plasma sheet during substorm phases observed by the interball/tail satellite     

M. Stepanova  T. Vucina-Parga  E. Antonova  I. Ovchinnikov  Y. Yermolaev   《Journal of Atmospheric and Solar》2005,67(17-18):1815

Fluctuations of the plasma bulk velocity across the plasma sheet are studied using single-point measurements from the Corall instrument on board the Interball/Tail satellite. Several hour-long intervals of continuous data corresponding to quiet geomagnetic conditions and different phases of isolated substorms are analyzed. The plasma sheet flow appears to be strongly turbulent, i.e. dominated by fluctuations that are unpredictable. Corresponding eddy diffusion coefficients were obtained as a function of the autocorrelation time and rms velocity of the fluctuations. It was found that the amplitude of the turbulence and the values of eddy-diffusion coefficients increase significantly during substorm growth and expansion phases and they decrease to their initial level during the recovery phase. We also studied a relationship between the eddy-diffusion coefficients and the absolute value of the geomagnetic field, also measured by the Interball/Tail satellite. It was found that this relationship varies depending on the phase of substorm, indicating possible change in the turbulence regimen with substorm phase.  相似文献   

Stochastic Modelling of the Impact of Flood Protection Measures Along the River Waal in the Netherlands     

Saskia Van Vuren  Huib J. De. Vriend  Sonja Ouwerkerk  Matthijs Kok 《Natural Hazards》2005,36(1-2):81-102

River flooding is a problem of international interest. In the past few years many countries suffered from severe floods. A large part of the Netherlands is below sea level and river levels. The Dutch flood defences along the river Rhine are designed for water levels with a probability of exceedance of 1/1250 per year. These water levels are computed with a hydrodynamic model using a deterministic bed level and a deterministic design discharge. Traditionally, the safety against flooding in the Netherlands is obtained by building and reinforcing dikes. Recently, a new policy was proposed to cope with increasing design discharges in the Rhine and Meuse rivers. This policy is known as the Room for the River (RfR) policy, in which a reduction of flood levels is achieved by measures creating space for the river, such as dike replacement, side channels and floodplain lowering. As compared with dike reinforcement, these measures may have a stronger impact on flow and sediment transport fields, probably leading to stronger morphological effects. As a result of the latter the flood conveyance capacity may decrease over time. An a priori judgement of safety against flooding on the basis of an increased conveyance capacity of the river can be quite misleading. Therefore, the determination of design water levels using a fixed-bed hydrodynamic model may not be justified and the use of a mobile-bed approach may be more appropriate. This problem is addressed in this paper, using a case study of the river Waal (one of the Rhine branches in the Netherlands). The morphological response of the river Waal to a flood protection measure (floodplain lowering in combination with summer levee removal) is analysed. The effect of this measure is subject to various sources of uncertainty. Monte Carlo simulations are applied to calculate the impact of uncertainties in the river discharge on the bed levels. The impact of the “uncertain” morphological response on design flood level predictions is analysed for three phenomena, viz. the impact of the spatial morphological variation over years, the impact of the seasonal morphological variation and the impact of the morphological variability around bifurcation points. The impact of seasonal morphological variations turns out to be negligible, but the other two phenomena appear to have each an appreciable impact (order of magnitude 0.05–0.1 m) on the computed design water levels. We have to note however, that other sources of uncertainty (e.g. uncertainty in hydraulic roughness predictor), which may be of influence, are not taken into consideration. In fact, the present investigation is limited to the sensitivity of the design water levels to uncertainties in the predicted bed level.  相似文献   

Improving the prediction of extreme FPSO hawser tension,using another highly correlated hawser tension with a longer time record     

《Applied Ocean Research》2019

The floating production storage and offloading unit (FPSO) is an offshore vessel that produces and stores crude oil prior to tanker transport. Robust prediction of extreme hawser tensions during the FPSO offloading operation is an important safety concern. Excessive hawser tension may occur during certain sea conditions, posing an operational risk. In this paper, the finite element method (FEM) software ANSYS AQWA has been employed to analyze vessel response due to hydrodynamic wave loads, acting on a specific FPSO vessel under actual sea conditions.In some practical situations, it would be useful to improve the accuracy of some statistical predictions based on a certain stochastic random process, given another synchronous highly correlated stochastic process that has been measured for a longer time, than the process of interest. In this paper, the issue of improving extreme value prediction has been addressed. In other words, an efficient transfer of information is necessary between two synchronous, highly correlated stochastic processes. Two such highly correlated FPSO hawser tension processes were simulated in order to test the efficiency of the proposed technique.  相似文献   

Retention of cosmogenic 3He in calcite     

《Quaternary Geochronology》2015

Cosmogenic ³He can be used to date a wide range of mineral phases because it is produced from all target elements and can be readily measured above atmospheric contamination. Calcite is a particularly attractive target mineral due to its natural abundance, large crystal size (>1 mm), and low He closure temperature (<70 °C), which limit non-cosmogenic ³He components (Copeland et al., 2007). However, several recent studies have shown that some calcite may not be retentive to helium, even under surface temperatures (Cros et al., 2014; Copeland et al., 2007). This study thus explores ³He retention and production in natural calcite samples at four different sites. Samples from two high elevation sites appear retentive to ³He over 10 kyr timescales, whereas two additional sites clearly suffer from diffusive loss of ³He. Step-degassing experiments suggest that diffusion in calcite is controlled by multiple diffusion domains, with an apparent activation energy of 25–27 kcal mol⁻¹. Although minor ³He loss is expected from the smallest diffusion domains, the observed kinetics cannot explain the poor retention at all sites. We thus propose that opaque (non-transparent) calcite may be more retentive due to the presence of imperfections in the crystal lattice. We conclude that ³He dating of calcite shows promise in some settings. However, because retention depends on crystallographic variability it must be evaluated on a case-by-case basis until robust criteria for retention can be identified.  相似文献